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ABSTRACT 

The dy nami cally zoned pulsation code developed by Castor, Davis, 
and Davison has been used to recalculate the Goddard model and to cal- 
culate three other Cepheid models with the same period (9.8 days). 

.This family of models shows how the bumps and other features of the 
light and velocity curves change as the mass is varied at constant 
period. This study, with a code that is capable of producing reliable 
light curves, shows again that the light and velocity curves 
for 9.8-day Cepheid models with standard homogeneous compositions do 
not show bumps like those that are observed unless the mass is signif- 
icantly lower than the "evolutionary mass." The light and velocity 
curves for the Goddard model presented here are similar to those com- 
puted independently by Fischel, Sparks, and Karp. They should be 
useful as standards for future investigators. 

Also Los Alamos Scientific Laboratory Consultant. 
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I. INTRODUCTION 

Cepheids with periods around ten days are known to show secon- 
dary bumps in their light curves around the time of maximum light 
(Cox 1974). Until recently, however, theoretical light curves for 
Cepheids have been unsatisfactory because of contamination by spurious 
"zoning bumps" (Keller and Mutschlecner 1971). The new dynamically 
zoned stellar pulsation code developed by Castor, Davis, and Davison 
(1977) is capable of producing reliable light curves without these 
artifacts. We have used a somewhat improved version of this code to 
recalculate the Goddard model and to calculate three other 9.8-day 
Cepheid models. We will discuss here the systematic trends that 
appear in the light and velocity curves for this family of models. 

II. THE MODELS 

The family of models represents a one-parameter sequence with 

mass as the independent variable, subject to the constraints that: 

1) the period is fixed at 9.8 days, and 2) the models lie along a 

line in the HR Diagram near the center of and roughly parallel to the 

instability strip. The characteristics of the models are given in 

Table 1. The actual masses of the models range from 62 to 99 percent 

of their respective "evolutionary masses" (corresponding to their 

luminosities with a standard mass-luminosity relation; see Table 1). 

% 

The light and velocity curves for the four models are shown in Figures 
1 and 2. 
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TABLE 1 


MODEL PARAMETERS 


Mass (M ) 
0 

4.0 

5.0 

6.0 

7.4 

Luminosity (L q ) 

3187. 

3849. 

4421. 

5210. 

Effective Temperature (K) 

5700 

5690 

5687 

5682 

& 

Evolutionary Mass (M q ) 

6.46 

6.82 

7.10 

7.44 

Zones in Model (Opt. Thin) 

65(12) 

72(15) 

73(12) 

74(10) 

Artificial Viscosity Coeff. + 

1.0 

2.0 

2.0 

2.0 

Fractional Core Radius 

0.12 

0.12 

0.13 

0.11 

Period (days) 

9.79 

9.78 

9.79 

9.78 

/ 9 

Peak KE (1(J ergs. Expanding) 

1.78 

3.30 

5.30 

9.70 

Periods Calculated 

46 

52 

37 

51 

Time Steps per Period 

1670 

1170 

1160 

1230 

Periods per Minute (CDC 7600) 

0.54 

0.67 

0.64 

0.67 

Comments 

Goddard 



Evolutionary 


Model 



Mass Model 


*Log [M(evolutionary)/M 0 ] = 0.287 log (L/L @ ) - 0.195. 
^Stellingwerf (1975) cutoff at 0.1. 

^Includes periods of amplification. 
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BOLOMETRIC MAGNITUDE 



Fig. 1. Bolometric light curves for four 9.8-day Cepheids. Zero phase 
is at maximum light, except for the 4.0 M 0 (Goddard) model, 
where zero phase was shifted 0.2 units earlier for easier com- 
parison. The model parameters are given in Table 1. 
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VELOCITY (KM/SEC) 
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Fig. 2. Velocity curves for four 9.8-day Cepheids. The velocities 

(positive for radial expansion) are interpolated values that 
show the motion of the matter at optical depth two— thirds, 
phase convention is as in Fig. 1. The model parameters are 
given in Table 1. 



The original version of the code, used by Castor, Davis, and 
Davison (in Deupree 1976; see also Davis and Davison 1978) to calcu- 
late their Goddard model, incorporated an exponential artificial 
viscosity cutoff. The principal change in the code since then has 
been to go to the standard Stellingwerf (1975) cutoff. . In addition, 
the value of the cutoff has been increased from one to ten per cent 
of the isothermal sound speed. These changes account for the higher 
limiting amplitude and the slightly different appearance of the light 
and velocity curves for the present Goddard model • compared with the 
earlier one (their model II). 

The light and velocity curves for the Goddar-d model in Figures 1 
and 2 should be useful as reference light and velocity curves for the 
Goddard model. They are quite similar to those presented by Fischel, 
Sparks, and Karp (in Deupree 1976). It is comforting that there is 
now reasonable agreement between the light and velocity curves calcu- 
lated by (at least two) independent investigators for the same model. . 

III. DISCUSSION OF THE LIGHT AND VELOCITY CURVES 

Let us first discuss the features in the light and velocity curves 
for the Goddard model. Starting at minimum light, we see the phase of 
rapidly rising light, interrupted by the so-called "artificial viscos- 
ity dip." The dip is associated with the compression wave that stops 
the infall of the envelope. The artificial viscosity plays a role in 
that it spreads out the compression wave in the mesh. However, the dip 
itself may not be entirely artificial, since the compression wave un- 
doubtedly raises the effective gravity briefly and can thus cause the 
light output to drop. 
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In this model, peak light is preceded by a shoulder in the light 
curve. There is a broad maximum in the velocity curve that begins at 
the phase of the shoulder and lasts through the light maximum. The 
light curve near maximum seems to be governed by changes in the effec- 
tive gravity at optical depth two-thirds, as there is an increase in 
the gravity between the shoulder and the peak and a gravity minimum at 
the peak itself. However, the physical origin of the shoulder is. not 
yet clear. 

After maximum light, there is a dip and then the "Christy bump,' 
which is associated with the reflection of a compression wave off the 
core. The bump in the velocity curve actually agrees in phase with 
the dip preceding the bump in the light curve. Therefore, in agreement 
with previous investigators, we identify the feature in the light curve 
as a "Christy dip," due to the compression that occurs when the re- 
flected wave reaches the surface. Although the light and velocity 
curves are similar, they differ in detail. This raises some questions 
about the practice of comparing observed light curves with calculated 
velocity curves. 

We now discuss how these features change as the mass is increased. 
In general, the Christy bump or dip moves down the descending branch. 

At the same time, it becomes less prominent, especially in the veloc- 
ity curve. In fact, there is no bump in the evolutionary mass (7.4 M q ) 
model. This confirms the conclusion that, with a standard homogeneous 
composition, the mass must be less than the evolutionary mass to show 
bumps like those that are observed. 
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Another interesting change that occurs as the mass is increased 
involves the light curve near maximum. In the 5.0 M q model, what was 
the shoulder in the Goddard model has become the peak, while the peak 
in the Goddard model has evolved into a bump following may-tn^ m light. 

This bump then moves down the descending branch and becomes less 
prominent as the mass is further increased. 

The ascending branch also changes with increasing mass. It be- 
comes steeper and the "artificial viscosity dip" moves up the ascend- 
ing branch until it nearly coincides with the peak in the evolutionary 
mass model. The compression wave associated with the dip increases 
greatly in strength as this happens. It should be emphasized, though, 
that this phenomenon is probably amplitude— dependent, and that the limit- 
ing amplitude, which also increases with mass (cf. Table 1), may depend 
on subtle physical or numerical effects in the calculation. 

Another interesting phenomenon occurred in the 4.0 and 5.0 M models 

© 

although it is not apparent in the light or velocity curves. In these 

models, at the phase of rapid expansion, when the hydrogen ionization 

front is racing inward, a shock forms on the neutral side following the 

front. This can be understood as a transition in the nature of the 

ionization front from D to R type (Kahn 1954, Castor 1966). Conditions 

change too rapidly for a steady-state R-type front to develop, and the 

front soon returns to its usual D— type behavior. This transition occurs 

about the time of the Christy dip in the Goddard model, and about the 

time of the dip just after maximum light in the 5.0 M model. The front 

© 

is near optical depth unity at this time, so there should not be any 
major changes in the emergent spectrum or energy distribution. However, 
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the following shock may be connected with the weak redshifted Ca II 
emission features that appear just after maximum light in Beta Doradus 
(Gratton 1953). The ionization front transition will be discussed in 
detail elsewhere (Adams and Castor, in preparation). 

IV. CONCLUSIONS AND FUTURE PLANS 

We have calculated a family of 9.8-day Cepheid models with a code 
that can produce reliable theoretical light curves. We find that only 
if the mass is significantly less than the "evolutionary mass" do the 
bumps in the light and velocity curves occur near maximum light, where 
they are observed to occur in real ten— day Cepheids (Cox 1974) • How- 
ever, only models with homogeneous population I compositions were cal- 
culated. It thus appears that we require either reduced masses, inhomo- 
geneous compositions, or some change in the physics in order to match 
the observed light curves of ten— day Cepheids. 

We plan to continue this program by calculating families of models 
at other periods to see how the light and velocity curves vary with 
period. Our goals will be to learn more about the mechanisms that pro- 
duce bumps in the light and velocity curves and, in general, to see if 
the Hertzsprung sequence (variation of bump location with period) can 
be reproduced theoretically. We would also like to calculate some in- 
homogeneous models to study Art Cox’s proposed solution to the "mass 
anomaly." 

Finally, we conclude that there is now a consensus standard Goddard 
model that can be used as a benchmark for any new pulsation code. It 
is entirely fitting that this consensus should be reported here at the 
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Goddard Space Flight Center only two conferences after the idea of a 
standard Goddard model was originated. 
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Discussion 


j. Wood : Do you have a temperature range for this 8 Dor Goddard model? From 

the G-band, I found about 1500° over the cycle, but Rogers and Bell found 
about 450° from Ha profiles. Sidney Parsons found about 900° from UBV work. 

Adams : First of all, I didn't make any attempt to match a particular star. . 

This was simply a theoretical sequence, so the parameters may not be correct 
for a particular star. The temperature data are available, and could be used 
to plot loops in the color-color diagram. However, there's a larger question 
that has to do with limiting amplitudes. The light amplitude shown here was 
rather large for a Cepheid; therefore, any behavior that depends on the 
limiting amplitude is questionable. But I don't know whether the temperature 
range depends on the limiting amplitude. 

Wesselink: Have you made any comparison between your light curves and observed 

light curves? 

Adams: I've just taken a cursory look at the observed light curves. I can't 

claim that I can match any observed light curve. The most important conclusion 
is that we can calculate good light curves, and they still show the traditional 
mass anomaly with the bumps. Let me add that even the features seen in the 
superb light curves produced by Pel and Lub still leave a little freedom for 
the theoretician to imagine various things happening. It's an extremely diffi- 
cult problem to get sufficient data. 

Connolly : What is causing that sharp dip at rising light and is it a dip, 

not actually a bump? 
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Adams: It is a dip. In the theoretical calculations it is due to the fact 

that when the envelope is falling in, it is stopped by a compression wave or 
shock wave. The intense compression for a short period of time steals some 
luminosity and puts it into ionization or compression. Therefore there is a 

sharp drop in the luminosity as the shock passes through the photosphere at 
x = 2/3. 

— sselink : I know one observed light curve — VZ Her — which resembles one 

of yours. It has a dip on the rising light. 

Adams: That is interesting. However, I would not go so far as to say that 
the dip has to exist theoretically, even though it shows up in the models. I 
think that it is difficult to find any observer who has a well-documented dip 
in the rising light because that's a very hard thing to see. It looks a lot 
like a cloud going by, lasting only a very short time. 

Pel: But they exist for longer period stars, with 15-30 day periods. There 

is a "standstill. 11 

Adams: That feature may be a bump, related to the bumps that define the 

Hertzsprung progression. 
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